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Abstract: During vitreoretinal surgery, the surgeon manipulates retinal 
tissue with tool-to-tissue interaction forces below the human sensory 
threshold. A force sensor (FS) integrated with conventional surgical tools 
may significantly improve the surgery outcome by providing tactile feedback 
to the surgeon. We designed and built a surgical tool integrated with a 
miniature FS with an outer diameter smaller than 1 mm for vitreoretinal 
surgery based on low-coherence Fabry-Perot (FP) interferometry. The force 
sensing elements are located at the tool tip which is in direct contact with 
tissue during surgery and the FP cavity length is interrogated by a fiber-optic 
common-path phase-sensitive optical coherence tomography (OCT) system. 
We have calibrated the FS's response to axial and lateral forces and 
conducted experiments to verify that our FS can simultaneously measure both 
axial and lateral force components. 
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1. Introduction 

Microsurgery involves the manipulation of very small anatomical structures. It is a series of 
manual manipulations performed predominantly using visual feedback, typically with the aid of 
a surgical microscope or other magnification device. In addition to visual feedback, tactile 
feedback can improve surgical outcome significantly. However, in microsurgical procedures 
such as vitreoretinal surgery, tool-to-tissue interaction forces are usually below human 
perception thresholds. Gupta et al. reported that 75% of tool-to-tissue interaction forces in 
vitreoretinal surgery were less than 7.5 mN in magnitude and that only 19% of force events at 
this magnitude were felt by the surgeon [1]. Although various tactile feedback or "sensory 
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substitution" schemes (e.g., [2-6]) may be used to provide feedback to the surgeon about these 
forces to help prevent potential tissue damage in freehand or robotic assisted surgery, a suitable 
force sensor (FS) must be incorporated into the surgical tools. In vitreoretinal surgery, the 
interaction forces between the tool and the sclera (white part) of the eye are often equal to or 
larger than the forces between the tool tip and the tissue being manipulated [7,8]. Thus, in 
addition to high sensitivity in force measurement, vitreoretinal microsurgery also requires using 
small size sensors that can be integrated into sub-millimetric instrument shafts. These 
considerations make the design of suitable force sensors for vitreoretinal microsurgery very 
challenging. 

Previously, a MEMS technique and miniature strain gauge have been implemented in the 
development of micro force sensing tools for retinal microsurgery, such as 1 degree of freedom 
(DOF) pick-like FS developed by Gupta et al. [1], multi-axis MEMS FS developed by Kim et 
al. [9], microgripper-integrated FS developed by Menciassi et al. [10], 3 DOF FS developed by 
Berkelman et al. [11,12] which had been incorporated later into a hand-held instrument and 
used to measure forces in retinal tasks [13]. Most of the above-mentioned sensors, with the 
exception of that of Gupta et al. [1], are too large to be incorporated within a surgical tool that 
will be placed inside the eye during vitreoretinal microsurgery. 

Force sensors based on fiber optics can be very small, mechanically stable, immune to 
electrical noise, and sterilizable [14]. Therefore they are ideal for incorporating into 
conventional surgical tools and for use as sensory substitution for tactile feedback in 
vitreoretinal microsurgery. Various force sensing schemes based on fiber optic techniques have 
been investigated over the last several decades. For example, FS that measures the change of 
optical power due to beam deflection has been developed for minimally invasive surgery 
[15,16]. However, suchFS cannot be directly used in vitreoretinal microsurgery due to the large 
diameter (about several to ten millimeters) and low force sensitivity (several to tens of 
millinewtons). In our previous work [8], we incorporated fiber Bragg grating (FBG) sensors 
into the shafts of microsurgical instruments. Although this approach produces very high force 
sensitivity in directions that are perpendicular to the tool shaft (lateral directions), the axial 
rigidity of both the tool and fibers makes achieving high axial force sensitivity with FBG 
sensors a challenging problem. Fabry-Perot interferometer (FPI) type sensors have also been 
used in bio-applications, including surgical force sensing [17,18]. The immediate advantage of 
FPI type micro force sensing is that it can potentially achieve a better axial sensitivity than that 
of FBG FS. 

In this study, we designed and built a fiber optic FS consisting of Fabry-Perot 
interferometer (FPI) cavities at the distal end of a fiber probe for potential applications in 
vitreoretinal surgery. In our FPI-FS, the length of FP cavity changed with force and we used a 
common-path Fourier domain optical coherence tomography (CP-FD-OCT) to interrogate the 
change of cavity length by phase-sensitive detection. The simple and compact common-path 
OCT configuration not only allows the use of optical probes with arbitrary lengths [19,20], but 
also provides improved phase stability to achieve sub-nanometer sensitivity in displacement 
measurement [21] and therefore to achieve high force sensitivity. Moreover, OCT also 
facilitates simultaneous measurement of multiple FP cavity lengths [22]. In our FS, we 
multiplexed three FPIs together to simultaneously obtain signals from three different channels. 
This enabled us to measure force in three dimensions with components both along and 
perpendicular to the tool shaft, which, to the best of our knowledge, has not been reported 
before. 

In this manuscript, we present the designing and fabrication of our FPI-FS (Section 2). 
Principle of phase-sensitive detection and multidimensional force measurement are presented 
in Section 3. In Section 4, we summarize the results of tool calibration and show experimental 
results of measured force. The final sections are discussion and conclusion. 
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2. Principle and tool design 

2. 1. System overview and tool fabrication 

Figure 1(a) shows the system schematic of FPI-FS for vitreoretinal microsurgery. The FP 
cavities that sense the forces are located at the tip of a conventional surgical tool — a surgical 
pick in this case — which is inserted into an eye to perform surgical tasks during vitreoretinal 
surgery. Detail of the FPI-FS probe is shown on the left side of Fig. 1(b). A CP-FD-OCT system 
is used to interrogate the length of the cavities by phase-sensitive detection. More details of our 
CP-FD-OCT system can be found in our previous work [23]. In brief, as shown on the right side 
of Fig. 1(b), it is a spectral-domain OCT system operating in 800 nm wavelength range with a 
~6um axial resolution and a ~3 mm imaging depth. The real-time 1-kHz data acquisition and 
processing was implemented with C++ in a Dell Precision T7400 computer. 

In our FPI-FS system, a superluminescent emission diode (SLED, Exalos Inc.) is used as 
broadband source (BS) to illuminate three FPIs which consist of a shared end mirror and three 
lead-in single-mode fibers (SM800-5. 6-125, Thorlabs, 125 urn cladding diameter). The mirror, 
which serves as the end reflector of the FPIs, was fabricated by polishing the surface of a 
stainless steel wire having 0.5mm diameter and 2.5mm length. We cleaved the fiber tips to right 
angle so the light is partially reflected at the fiber tip. Optical fields reflected by the fiber tips 
(E u ,i = 1,2,3) and the metal reflector (E 2 -„i = 1,2,3) are combined (as shown in Fig. 1(b)), routed 
through the coupler, and the interference signal is detected by the spectrometer. 

As depicted in Fig. 1(b), the force-sensing part of the tool contained three optical fibers and 
a mirror, which form three FPIs with different cavity lengths. Other components of the FS 
include a tool shaft, a flexure, and a tool tip (pick). The tool shaft is made of a piece of titanium 
wire that is strong, lightweight, and biocompatible. The wire is 0.5 mm in diameter and 50 mm 
in length. Three square section channels (0.15 mm x 0.15 mm) were machined (laser cut) into 
the surface of the wire along the tool axis, 2jt/3 from one another in a circle, as shown in Fig. 
1(c). Fibers were embedded and fixed into the channels. 

The mirror and tool shaft are bound by a flexure made of Nitinol using a laser cut technique. 
The flexure deforms proportionally to the force applied and therefore induces change in the 
lengths of the FP cavities that will be used to calculate the force. We chose to use Nitinol 
because of its super-elastic property which enables the tool to have a measurable deformation 
with a force in the order of a millinewton. Nitinol has modulus of elasticity ranging from 41 x 
10 3 MPa to 75 x 10 3 MPa; therefore, the fabricated tool has extremely high force sensitivity. As 
shown in Fig. 1(d), the flexure — with 0.80mm outer diameter (OD) and 0.60mm inner diameter 
(ID) — has five layers; each layer consists of three 40-micrometer struts; and between adjacent 
layers there is a 100 micrometer slot. The multi-layered configuration can further increase the 
elasticity of the flexure in the axial direction and therefore achieve a high sensitivity for force 
measurement. The total length of the flexure is 3 mm, with 1 mm between the distal end and the 
first slot. The pick was made using a stainless steel SS304 hypodermic tube, 21 gauge (OD = 
0.80mm, ID = 0.55mm). The tip was bevel cut (at 75 degrees) and bent at around 45 degrees 
relative to the axis. The pick total length is 2.5mm, with a 1.5mm length bevel section. Loctite 
medical instant adhesive (Henkel, Dusseldorf, Germany) was used to bond all the parts 
together: first the mirror with the pick, next the flexure, and finally the tool shaft. The optical 
fibers were embedded into slots machined to the tool shaft. The distance between the fiber tips 
and the polished metal surface was manually adjusted before finally bonding the fibers with the 
tool shaft. It is worth mentioning that we made the FP cavities to have different lengths. This 
allows us to detect the length change of each cavity respectively. The CAD model for the 
FPI-FS probe integrated with the pick is shown in Fig. 1(e). For calibration and future manual or 
robotic use, the force sensing probe is attached to a custom-made handle, as shown in Fig. 1 (f). 
Figure 1(g) is a photo of the tool pictured below a U.S. quarter. 
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Fig. 1. (a) System schematic of FPI FS; (b) detailed schematic of the system including FPI-FS 
and CP OCT interrogation; (c) cross section of tool shaft with fiber embedded; (d) CAD model 
for the Nitinol flexure; (e) CAD model of the FPI-FS probe without the tool handle; (f) CAD 
model of the FPI-FS probe with the tool handle; (g) photo of the tool pictured below a US 
quarter. 

2.2. OCT signal and phase-sensitive detection 

Assuming that the polished metal surface reflects all the incident light and that the fiber tip 
reflectivity, r, is about 4% according to Fresnel equation [24], we calculated the finesse of the 
FPI to be about 1.25. The round trip loss of optical power in the FP cavity is larger than 4% due 
to beam divergence; therefore, the actual finesse of FPIs are even lower than 1.25. Moreover, 
we have measured the fringe visibility of interferograms obtained from the FPIs and the values 
are very small (0.39, 0.047 and 0.049, for three FPIs, respectively). Therefore, in our 
multiplexed FPIs, signals due to higher order reflections are small and therefore can be 
neglected in the following analysis. 

As shown in Fig. 1(b), we denote the optical field reflected by the tip of the i th fiber as E 1; . 
We denote the optical field reflected by the shared end mirror and collected by the i fiber as E 2 , 
(z = 1,2, and 3). We have E 1; = rajE 0 (k)exp(-jklj) and E 2 , = r 2 ia i E 0 (k)e\p[-jk(l i + 2Lj)]. Here,y" is 
the imaginary unit; k = 2n/X is the wavenumber; E 0 (k) is the spectral density of the broadband 
source; a, is the portion of light that is coupled into the i th fiber from the source; r 2 , takes into 
account both fiber coupling loss and the reflectivity of the polished metal surface; /, is the round 
trip optical path length traveled within the ; th fiber; Z, is the length of the i FP cavity, which 
varies when there is force applied to the tool, as shown in Eq. (1): 

L, =L uo +Sl i (F). (1) 
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In Eq. (1), L i 0 is the length of i* FP cavity at neutral when there is no force exerted. 5/, indicates 
change of cavity length due to force exerted to the tool. 

Since the round trip optical path length traveled within the / th fiber, /,, varies significantly 
from other fibers (much more than several millimeters), Ey or E 2i does not interfere if vtj. Based 
on this fact, the spectral interferogram S (k) can be expressed as in Eq. (2): 



S(k) = 77 



^\E u (k)\\^\E 2i (kf + ^2Rc(E u (k)E 2 ;(k)) 



(2) 



In Eq. (2), rj indicates the responsive coefficient of the system; * denotes complex 
conjugate; Re() refers to the real part of a complex signal. The third term of the right-hand side 
of Eq. (2) is the interference term. OCT signal, I OC t(z), can thus be obtained by performing 
inverse Fourier transform on S(k): 



I OCT (z) = F- l [s(k)] (3) 



Denoting the axial point spread function of the OCT system as h(z), which equals 
F~ l (\E Q (k)^\ we can rewrite Eq. (3) as Eq. (4) with approximation when the bandwidth is 
relatively small: 



I oct( z )= 7M z )ZfKT + K a < 

+ n f j [r r ; i afh(z-2L. 0 )e j2 ^] (4) 

i=l 

3 

+vY\_ r ' r 2< a < h { z+2L <*) e ~ 



-jlk 0 Sl, 



The first term on the right-hand-side (RHS) of Eq. (4) is the auto-correlation term. The second 
and third terms in the RHS of Eq. (4) represent coherence peaks due to the interference between 
two end reflectors of each FPL The second and third term in the RHS of Eq. (4) are symmetry 
about the zero delay because the Fourier transform is performed on real spectral data. 

As indicated by Eq. (4), the change of FPI cavity length 5/, can be extracted from (j) h the 
phase of OCT signal at delay L ifi , which equals tan" 1 {Im[/ OC7 {Z, ii0 )]/ Re[I 0C7 (L ifi )]}. However, 
the function tan -1 gives value in the range of [-n/2, rc/2]; therefore 

2^. = tan--{ 4^4 }-^ (5) 

4* [Re[/ OCT (A,o)]J '4 

In Eqs. (5) and (6), N t is the integer that makes (5/, + NM4) to fall within the range of [- V8, 
A.o/8]. With small force variation, A^- stays constant. Denoting the first term in the RHS of Eq. (6) 
as di and the second term as d i?0 which is assumed to be a constant, we have 

Sl i (F) = d i -d. 0 (7) 
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When the variation of FP cavity length is large — such that $ jumps from -nil directly to nil or 
from ti/2 directly to -nil — a phase unwrapping technique can be used to obtain continuous 
phase or displacement [21]. 

When a small force is applied to the tool, the change in FP cavity length is extremely small 
(in the order of a nanometer) compared to the axial resolution of our OCT system — even when 
the super-elastic material Nitinol is used to fabricate the flexure. Therefore it cannot induce any 
observable change in the amplitude of OCT signal. Nevertheless, the small displacements can 
lead to a measurable phase shift of the interference spectrum. In other words, the small 
deformation can be extracted from the phase of OCT signal at the coherence peaks as shown in 
Eq. (6), which is demonstrated in the following simulation. In this simulation, we assumed Zi 0 
= lOOurn, Z 2 ,o = 150 urn, and Z 3 0 = 200um; the light source has central wavelength of 800nm 
and FWHM bandwidth of 60nm. Superposed spectral interference signals from three FP 
cavities are simulated and shown in Fig. 2(a) according to Eq. (2). Performing inverse Fourier 
transform on this interference signal leads to complex-valued spatial domain OCT signal 
^oct(z), which contains three coherence peaks at z = L u L 2 , and L 3 corresponding to the length of 
the three FP cavities. The amplitude of OCT signal is shown in Fig. 2(b). With force exerted, the 




wavelength (nm) delay 



Fig. 2. (a) Superposition of spectral interference signals from three FP cavities; (b) OCT signal 
with three coherence peaks corresponding to three FP cavities; (c) central part of the interference 
spectrum: with (red) and without (black) additional displacement <5/,; (d) complex OCT signals at 
the 1 st, 2nd, and 3rd coherence peaks in the complex plane using polar coordinate system. Black 
symbols and red symbols represent OCT signals with and without additional displacement <5/;; (e) 
actual interferometric spectrum obtained from the FPI-FS; (f) amplitude of OCT signal obtained 
by Fourier transforming spectrum shown in Fig. 2(e). 
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cavity length changed proportionally to the force. In this simulation, we assumed 8l\ = 20nm, 
8l 2 = 40nm and c*/ 3 = 60nm. However, such small change in cavity length cannot be seen directly 
from the amplitude of OCT signal described in Fig. 2(b). On the other hand, as shown in Fig. 
2(c) which is the central part of the interference spectrum, change in cavity length in the order of 
nanometer shifts the original spectrum (black) to the red one. To demonstrate that the small 
displacement can induce measurable phase change to the complex-valued OCT signal at the 
coherence peak, we show complex OCT signals at the 1st, 2nd, and 3rd coherence peaks in the 
complex plane using a polar coordinate system in Fig. 2(d). Black symbols and red symbols 
represent OCT signals with and without additional displacement dk. Clearly, the vector that 
represents the complex OCT signal rotates due to phase shift induced by 81,. 

In Fig. 2, we also show the actual spectrum obtained from the multiplexed FPIs of our FS as 
Fig. 2(e). We performed Fourier transform on the interferogram and obtained OCT signal 
shown in Fig. 2(f). Three peaks in Fig. 2(f) come from interference between end reflectors of 
the three FP cavities and the locations of these peaks indicate the length of FP cavities. 

2.3. 3D force measurement 

As demonstrated in the three-dimensional Cartesian coordinate of Fig. 1(b), a force (F) along an 
arbitrary direction can be decomposed into axial and lateral components. In the force coordinate 
attached to the tool tip, z axis or axial direction is along the tool shaft. We denote force along the 
tool shaft direction as axial force F z . x-y plane is perpendicular to z axis and force in x-y plane is 
denoted as lateral force F/. The choice of x andj axes is arbitrary as long as x, y and z axes form 
a right-hand coordinate system. 

As shown in Figs. 3(a) and 3(b), with an axial force applied to the tool, the flexure deforms 
and the change of FP cavity length is denoted as (5/ i z which equals (F Z D)/(A 0 E) [25]. Here D is 
the length of cantilever beam; E is the effective Young's modulus of the cantilever beam; A 0 is 
the cross-sectional area on which the force is exerted. When lateral force is applied, the beam 
deflects with an angle a that equals (FiD 2 )/(2EI) where I is the area moment of inertia of the 
wire cross section [25]. We show the beam deflection with exaggeration in Fig. 3(c). In fact, a is 
extremely small due to the small force exerted during vitreoretinal microsurgery. Therefore, the 
slight change in light beam propagation direction does not reduce the coupling efficiency when 
light gets reflected back to the single-mode fiber and does not reduce the OCT signal amplitude 
either. 

According to Fig. 3(d), the change of FP cavity length dkj due to lateral force can be 
expressed as [-2Z,, 0 sin 2 (a/2)/cos(a) + 6,tan(a)]. As a is extremely small, the following 
approximation is valid: sin(a/2)«a/2; tan(a)«a; cos(a)«l. Moreover, sin 2 (a/2) is a higher order 
term of the already extremely small value a and therefore can be negligible compared to tan(a). 
As a result, Shj = Z>,a. Here, Z>, is the distance between the center of the i th fiber and the neutral 
surface which is perpendicular to the applied force and passes through the center of the tool 
shaft cross-section. Denoting the angle between vector F/ and x axis as 9 as in Fig. 3(e), we find 
that bi equals (x/COsG-j/sinG) where (x,,y,) is the transversal coordinate of the i th fiber. Therefore 
bi can also be expressed as: 



In Eq. (8), e, = tan" (y/x;). Assume that the cavity change induced by force F is the linear 
superposition of the effect of F z and F t . In other words, the cavity length change of the z* FP 
cavity 8l\ equals c*/ i z + dly and we can express 8l\ as follows: 




(8) 
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F : D 



F +x. F, cos&-y F,sm& 

\E 2EI 2EI 



(9) 



It is clear in Eq. (9) that the displacement induced by lateral force has a sinusoidal 
dependency on 9. 





i 

























Fig. 3 . (a) FS at neutral; (b) axial force changes the length of flexure and changes the FP cavity 
length; (c) lateral force leads flexure to bend and changes the FP cavity length; (d) diagram 
shows the calculation of cavity length induced by lateral force; (e) diagram shows the calculation 
ofb, 

Denoting F/sinO and 7 7 /COs9 as F x and F v ; D 2 Xj /(2EI) and £) 2 j, /(2EI) as A ix and A iy ; DI{AqE) as 
A iz , we can re-write Eq. (9) as 



SI = A /•; + AF + A, / • 



(10) 



After introducing F x and F r , SI, is linearly dependent on F x and F y with constant coefficients 
A ix and^4„, which are independent of the direction of applied lateral force, as shown in Eq. (9). 
Therefore, A ix , A jy and A j: can be obtained from a linear regression procedure if we know the 
force applied to the tool and the corresponding displacements. 

Although for a given amplitude of lateral force F h F x and F y take different values if 9 varies, 
we can always calculate F t with Eq. (1 1), which is independent of 9: 



f,=JfJ7f} (ii) 

With Eqs. (7) and (10), the relationship between 61 \ and F can be shown as follows: 

(12) 
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Values of elements in matrix^ can be obtained from a calibration procedure shown in detail 
in the following section. Once A is obtained, we are able to calculate the force in three 
dimension using Eq. (13): 



F = A-'Sl 



(13) 
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As shown in Eq. (12) and Eq. (13), it is critical to subtract d 0i from the displacement c/, to 
obtain correct force measurement. Above do, is the displacement measured from the phase of 
OCT signal when no force is applied to the FS. If the subtraction or biasing procedure is taken 
with a certain force Fq applied to the FS, the measurement provides force relative to Fq rather 
than absolute force. 

3. Results 

3.1. Axial calibration 

In our axial calibration experiment as illustrated in Figs. 4(a) and 4(b), we attached different 
numbers of testing weights to the FS and maintained the tool shaft in the gravity direction to 
keep F x and F y equal to 0. Therefore, pure axial force along the tool shaft was applied in this 
setup and coefficients A iz could be obtained. In our calibration experiments, each testing weight 
was 0.897 ( ± 0.005) mN. 

Using the complex-valued OCT signal acquired from our CP-FD-OCT, we calculated the 
displacements at the three coherence peaks from the phase of OCT signal with Eq. (6) using TV, 
= 0. The measured displacements are equal to change in the length of FP cavities. With different 
numbers of testing weights i.e. axial forces, we obtained d\, d 2 , and <i 3 , displacements 
corresponding to the three FP cavities, shown as black circles in Figs. 4(c), 4(d) and 4(e), 
respectively. Afterwards, with the known axial force F : and the measured displacement d, = d 0i 
+ A iz F z , we were able to perform linear regression to extract A lz , A 2z , and A lz . The linear fitting 
results are shown as red lines in Figs. 4(b), 4(c) and 4(d). 




force (mN) force (mN) force (mN) 



Fig. 4. (a) Photo of axial calibration setup; (b) schematic of axial calibration setup; (c) 
displacement versus axial force (first FP cavity); (d) displacement versus axial force (second FP 
cavity); (e) displacement versus axial force (third FP cavity). 

3.2. Lateral calibration 

In the lateral calibration setup as shown in Figs. 5(a) and 5(b), we maintained the tool shaft (z 
axis of the force coordinate) perpendicular to the gravity direction so that the force induced by 
the gravity of the testing weights had only a lateral component in the force coordinate that was 
attached to the tool tip; in other words, F z = 0. As indicated by Eq. (9), the displacement induced 
by lateral force not only depends on the magnitude of lateral force exerted to the FS, but also 
depends on 8, which essentially is the direction of lateral force in x-y plane as shown in Fig. 3 
(e). Therefore, to characterize the tool's response to lateral force, we needed to apply force in the 
lateral plane {x-y plane) with different magnitudes as well as at different azimuth angles 8. In 
order to do this, we attached the tool to a rotary stage in x-y plane to change the angle between 
the x axis of our force coordinate and the direction of gravity. Since our force coordinate was 
attached to the tool tip, 8 changes in the same manner as the rotation of the tool. 
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In our lateral calibration experiments, we attached testing weights to the tool; we rotated the 
tool circumferentially to obtain different values of 9 and took the displacement measurements. 
We repeated the above procedures with different lateral forces by attaching different numbers 
of testing weights to the FS. The experimental results of our lateral calibration are shown in 
Figs. 5(c)-5(e). 

Displacements from the 1st, 2nd, and 3rd FP cavities are indicated by circles in Figs. 
5(c)-5(e), in which data points with different colors were obtained with different lateral forces 
(numbers of testing weights are shown as legends in these figures). Figures 5(c)-5(e) clearly 
indicate that the displacements induced by lateral force vary sinusoidally as the tool orientation 
0 for a given lateral force, which is consistent with Eq. (9). Sinusoidal fittings of the 
experimental results are also shown as dashed curves in Figs. 5(c)-5(e) to further verify the 
functional dependency of displacement on 9. In addition, results in Figs. 5(c)-5(e) also show 
that the phase of the displacement's sinusoidal dependency on 9 varies for different FP cavities. 
This is because the lead-in fiber for the z* FP cavity takes a different transversal coordinate 
(x,,j,) and therefore the phase e, = tan _1 (x,/y,). Moreover, our results also show that the 
displacements increase proportionally as lateral force increases. To demonstrate this more 
clearly, in Figs. 5(f)-5(h), we show displacements obtained with different lateral forces when 9 
was 3n/2. A linear relationship between displacement and lateral force can be observed in Fig. 
5(f)-5(h). 

To obtain A ix and A- K using data shown in Fig. 5(c)-5(e), we performed a two-dimensional 
linear regression to solve the linear model d = A/F xy . Here A/ = [A hx A lv ; A 2x A 2y ; A ix A^]. F xy 
indicates known vectors with F x and F r as elements: F xy = [F X ; F y ]; F x = F/ cos9; F y = F/ sin9; F/ 




Fig. 5. (a) Photo of lateral calibration setup; (b) schematic of lateral calibration setup; (c)-(h) 
results from lateral calibration: (c)-(e) show displacements measured from three FP cavities at 
different 9 with different lateral forces. Legends indicate number of testing weights applied to the 
sensor; (f)-(h) show displacements with different lateral loads when 8 was 3jr/2. 
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is the known gravity of testing weights, d indicates the vectors with measured displacements 
from different FP cavities as elements d = \d\, d 2 ; c/3]. 

With At obtained from linear regression of the lateral calibration data, we could calculate F/ 
inversely. F x and F y were obtained by solving the linear equation d ~ A t F xy in a least-square 

t- IT 1 7 

manner: F xy = (A/ Ai) A\ d. Here () indicates to take the inverse of a matrix and () indicates 
to take the transpose of a matrix. F x , F v , and F/ extracted are shown in Fig. 6. Only two curves 
are visible in results with 9 = 0 degree, because F y is almost 0 and therefore F\ is almost 
overlapped with F x . Similarly, only two curves are visible in results with 9 = 90 degree, because 
F x almost is 0. Clearly, although F x and F v depend on 9 for the same lateral force, F calculated 
from Eq. (13) is independent of tool orientation, shown as black curves in Fig. 6. 



0 degree 



45 degree 



1 35 degree 




applied F (mN) 



applied F (mN) 



applied F (mN) 



applied F (mN) 



Fig. 6. F„ F r and F t extracted from displacements when 9 took different values. 

In conclusion, we obtained the following coefficients as elements of our calibration matrix 
A:A U = -30.3nm/mN; A ly = 27.5nm/mN; A Xz = 4.3nm/mN; A 2x = 33.9nm/mN; A ly = 8.9nm/mN; 
A lz = 6.2 nm/mN; A 3x = -1 1.8nm/mN; A 3y = -30.9nm/mN; A 3: = 6.5nm/mN. A ix = D 2 x t /(2EI) 
and Ajy = D yt /(2ET); and x, and y-, can take different values, positive and negative as shown in 
Fig. 3(e); therefore, A ix and A iy can take different value and have different signs for the same 
channel. Force could thus be calculated using Eq. (13) with A and the measured displacements. 

3.3. 3D Force measurement 

To validate that our FS can measure force with both axial and lateral components, we conducted 
the following experiment. In addition to the rotary stage in x-y plane, another rotary stage was 
used to change tp, the intersection angle between the tool shaft and the gravity direction, as 
shown in Figs. 7(a) and 7(b). With the same load attached to the FS, different values of <p leads 
to different values of axial and lateral force. Here G is the gravity of the load which equals 
8.2mN for this experiment. We used a rotation stage (as indicated by Rotation 1 in Figs. 7(a) 
and 7(b)) to change <p and therefore change F x , F y , F z correspondingly, because: 



F x = Gs'm<pcos8;F v = Gs'm(ps'm8;F : = Gcosip (14) 

We repeated the measurements with 9 = n/2 and 9 = Jt. Axial and lateral forces obtained are 
shown as circles in Figs. 7(c) and 7(d), respectively. Red circles indicate data obtained with 9 = 
7T./2; blue circles indicate data obtained with 9 = n. 
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It is worth mentioning that we biased the displacement reading when tp = n/2 and the FS was 
exerted with a non-zero force. As shown in Eq. (7) and discussed at the end of Section 2.3, the 
forces measured were values relative to the force when biasing was taken. As a result, the 
measured forces are F xm = F X {(p)-F x (<p% = M ; F yM = F^cp)-F^(p)\ v = na ; F z<m = F \((p)-F \{<p% = M - 
Incorporating Eq. (14), we have F/,,, = G (1-sinp) and F :m = Gcos<p, which are also shown in 
Figs. 7(c) and 7(d) as black curves. The consistency between force measured from our FS and 
the calculated black curves implies that our FS is able to measure force with both axial and 
lateral components. With the known load gravity (G = 8.2mN), we may bias the data in post 
processing and obtained the results shown in Figs. 7(e) and 7(f). Results in Figs. 7(e) and 7(f) 
show maximum axial force and minimum lateral force with tp — 0; as well as minimum axial 
force and maximum lateral force with tp = n/2. 




angle (degree 



angle (degree) 



Fig. 7. (a) Photo of 3D force measurement setup; (b) Schematic of 3D force measurement setup; 
(c) Measured and actual axial force; (d) Measured and actual lateral force; (e) Measured and 
actual axial force obtained by using the known load gravity as bias; (f) Measured and actual 
lateral force obtained by using the known load gravity as bias. 

4. Discussion 



Measuring a force in an arbitrary direction with respect to the tool shaft of the FS requires at 
least three FP cavities to provide independent displacement readings. First, the change of FPI 
cavity length is a linear superposition effect of axial and lateral forces. As indicated by our 
calibration experiments, the same magnitude of lateral force can induce much larger change in 
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cavity length compared to axial force. In other words, our FS's response to axial and lateral 
forces are very different. Therefore, it is necessary to differentiate F z and F/. Second, the 
displacement introduced by lateral force F t is angle (8) dependent, as shown in Eq. (9). As a 
result, it is more convenient to decomposition F/ in x-y plane as F„ F y and afterwards use Eq. 
(11) to calculate F h which would be independent of the choice of x or y axis of our force 
coordinate. As we need to determine three unknowns — F„ F y , and F z — it requires at least three 
independent linear equations based on measurements from three different FP cavities. 

As discussed in Section 2.3, we can directly calculate our FS's response to axial and lateral 
forces using the mechanical and geometric properties of our tool. However, the mechanical 
property of flexure might change significantly because Nitinol changes its mechanical property 
during laser cutting procedures. Moreover, we applied medical adhesive to glue different 
components together, which would also change the overall mechanical properties of the tool. 
Besides, due to the extremely small dimension of the tool, the geometry of the fabricated tool 
might be different from the initial design. As a result, experimental calibration can characterize 
the performance of our FS much better than a theoretical model or even finite element analysis. 

Our calibration results show that different FP cavities have different responses to axial force 
(A iz = 4.3nm/mN; A 2z = 6.2nm/mN; A 3z = 6.5nm/mN). However, as indicated by Eq. (9), the 
sensitivity of our FS to axial force should be identical for different FP cavities because A iz = 
DI(A<JE). The difference in axial sensitivity might be due to the asymmetry in mechanical 
property of our FS. Similarly, lateral sensitivity (A ix 2 + A iy 2 ) m = D 2 r/(2EI) should be the same 
for different cavities as in Eq. (9); however, our experimental values are 43.2nm/mN, 
35.0nm/mN, and 33.1nm/mN for the first, second, and third FP cavities, respectively. The 
difference in the lateral force sensitivities might be the result of the asymmetry in mechanical 
property of our FS, as well as the misalignment of the fibers. Due to the extremely small scale of 
our tool, it is almost impossible to make sure that the libers have the same distance to the center 
of the tool shaft to make (A ix 2 + A iy 2 ) 112 identical for the three FPIs. 

The force measurement sensitivity depends on the sensitivity of phase-resolved 
displacement measurement as well as A, coefficients relating cavity length change, and force. 
According to Eq. (13), the variance of force a F 2 = Var(F) can be calculated from the covariance 
matrix of Gaussian random vector SI, co\(dl): co\(F) = ^cov^OK^ -1 ) 7 ], and a F is the 
diagonal elements of the matrix cov(F) [26]. To estimate the magnitude of a F , we further 
assume 8l\, 5l 2 , <5/ 3 are independent random Gaussian variables with the same variance of. This 

2 2 2 1 IT 

simplifies the expression of a F , which becomes a F =a t A [(A ) ]. Although the fundamental 
lower limit of a t is determined by the signal-to-noise ratio (SNR) of the OCT system [17], phase 
or displacement fluctuation due to environmental disturbances such as temperature variation 
and air flow is slower but larger in magnitude, and therefore determines the sensitivity of force 
measurement. In our experiment, 07 was 0.5nm from a long-term (10s) measurement while 07 
was 0. lnm from a short-term (0.5s) measurement. Using 07 = 0.5nm and the obtained calibrating 
matrix A, we were able to calculate cov(F) and further extract a F : er^O.OlmN; <r fjI «0.01mN; 
ff ft «0.05mN. There results indicate that our FS can achieve sub-millinewton sensitivity in force 
measurement. 

As shown in Eq. (5) and Eq. (6), we extract displacement from the phase of a complex 
number that is defined in the domain [-jt/2 jt/2]. When a time -varying phase reaches and crosses 
the boundary of this range, the phase calculated from Eqs. (5) and (6) will show an abrupt 
change of ± n, depending on whether the phase decreases or increases. According to Eq. (7), the 
value of di is from -lOOnm to +100nm and therefore we can obtain displacement without 
ambiguity only in a range of 200nm. Since our FS has 30 -40 nm/mN sensitivity to lateral force, 
a 200nm displacement range allows us to measure 5mN ~6mN lateral force at most. Therefore, 
a severe restriction in the dynamic range of our force measurement is placed by the limited 
range of displacement measurement. However, we implemented a well-established method 
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called phase unwrapping in this work to increase the dynamic range of our force measurement 
[16]. In our C++ software, we continuously calculate the phases corresponding to each 
coherence peak and check whether there is an abrupt change of ±jt between two consecutive 
phases: (f) w <j) n +\. If such phase discontinuity occurs, we modify the phases obtained after <f> n by 
subtracting or adding %. This method can increase the dynamic range in the measurement of 
displacement and thus force, because the force is assumed to be continuous and usually cannot 
lead to such a big change in <f> due to the short time between two sampling points (1ms). 
However, if the force varies quickly (JF/dt>5mN/ms), the unwrapping technique cannot 
provide correct force measurement. 

5. Conclusion 

In this study, we proposed and studied a miniature fiber-optic force sensor based on low 
coherence Fabry-Perot interferometry for vitreoretinal surgery. We used a common-path 
Fourier domain optical coherence tomography to interrogate the change of FP cavity length 
through phase-sensitive detection. We fabricated a force sensing tool with an outer diameter 
less than 1 mm. Calibration was performed to characterize the tool. We also conducted 
experiments to verify that our FS can measure force with both axial and lateral components. 
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